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Abstract. We give a new explicitly invertible approximation of the normal cu- 
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mulative distribution function: $(x) ~ \ + \ \J 1 - e x 26.694+2*2 , Vx > 0, with 
absolute error < 4.00 • 10~ 5 , absolute value of the relative error < 4.53 • 10~ 5 , 
which, beeing designed essentially for practical use, is much simpler than a pre- 
viously published formula and, though less precise, still reaches 4 decimals of 
precision, and has a complexity essentially comparable with that of the approxi- 
mation of the normal cumulative distribution function $(2;) immediatly derived 
from Winitzki's approximation of erf (2), reducing about 36% the absolute error 
and about 28% the relative error with respect to that, overcoming the threshold 
of 4 decimals of precision. 
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This paper is devoted to approximate some special functions, in particular the 
normal cumulative distribution function. 

Though computers now allow to compute them with arbitrary precision, such 
approximations are still worth for several reasons, including to catch the soul 
of the considered functions, allowing to understand at a glance their behaviour. 
Let's add that, despite technologic progress, those functions - of wide practical 
use - not always are available on pocket calculators. In practice, the ancient 
numerical tables are still widely used; but they give approximations only for 
some values and, if linearly interpolating to approximate intermediate values, 
both precision and simplicity of use reduce. And surely by computers you may 
obtain graphs of those functions, but for a mathematician the meaning content 
of formuls is greater. 

Furthermore, here we produce an explicitly invertible (and, in fact, simply) ap- 
proximation, which allow to keep coherence working contemporarily with the 
considered function and its inverse. 



For the [5] [5] special functions 

^=e^"dt, (1) 
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and the related error function erf(x) — 2$(a;v / 2) — 1 (Vx € M), there are sev- 
eral approximations; in particular see classical [1] [3] [4] and recent [2] [5] [7] 
[5]; approximations for x > are sufficient because of the symmetry formula 
$(-sg) = 1 - $(x) (Vz e JR). 

Restricting now our attention only to those approximations which are simply 
explicitly invertible - in the sense, explicitly invertible without requiring to solve 
cubic or quartic equations - the most precise appears to be [5] 



1 1 [ -1.2735457*2-0.0743968x4 ( \ £ {x)\ < 1.14 . lO" 5 

$ ( , )S _ + _Vl-e~ !.2 + ,_4 | |^ 1<L78 . 10 -5 Vx>0 

(2) 

which is an improvement preserving (despite the adding of the quartic mono- 
mial) the simple explicit invertibility (essential solving a biquadratic equation 
after obvious substitutions) of this approximation of $ 
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\ 1 l \L — , 2/ f \e(x )\< 6.21 • 1(T 5 w ^ n , , 

$(a;)~ + Vl-e 2 (1+0 0735 ^ LV; L fi on.in-5 > (3) 



2 2 V \ \e r (x)\ < 6.30 ■ 10 

immediately derived by = ^ + 5 er /(^7|) (Va: > 0) from this [7 Winitzki's 
Approximation of erf 



r - t , -^i±°^ f |e(a;)| < 1.25 • HT 4 w ^ n 
er/(«) a Y 1 - e — | j^, < L28 . 10 - 4 Vx > . (4) 

In this note we give this new (simply) explicitly invertible approxi- 
mation of the normal cumulative distribution function 



*/ n i , ufi 17+ * 2 2 f |£(x)| < 4.00- 10" 5 w 

which, beeing designed essentially for practical use, 

• is much simpler than ([2]) and, though less precise, still reaches 4 
decimals of precision; 

• has a complexity essentially comparable with that of ([3]) reducing 
about 36% the absolute error and about 28% the relative error with 
respect to that, overcoming the threshold of 4 decimals of precision. 
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Instead, the corresponding approximation of erf is not so worth, because, though 
reduces about 36% the absolute error of Q, it remains with the precision of 
3 decimals, and furthermore the absolute value of the relative error |ey(o;)| < 
1.79 • 10~ 4 is quite greater than in Q. 

See below the graphs (made by Mathematica ®) of the approximation, of the 
absolute error and of the absolute value of the relative error for < x < 7. 
For x > 7 the trivial approximation $(a;) ~ 1 has absolute error and absolute 
value of the relative error highly less than 4 ■ 10 -5 and 4.53 • 10~ 5 respectively. 
(Nevertheless, if interested in a formal proof of the majorization \s(x)\ < 4-10 -5 
for x > 7, you may follow [S]). 
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Fig. 1. The new approximation ^ 
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Fig. 4. Second zoom of Fig. 2. 
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Fig. 2. Absolute error Fig. 5. Absolute value of relative error 
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Fig. 3. First zoom of Fig. 2 
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Fig. 6. Zoom of Fig. 5. 
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